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Abstract 



We discuss the optimization of optical microcavity de- 
signs based on 2-d photonic crystals for the purpose of 
strong coupling between the cavity field and a single neu- 
tral atom trapped within a hole. We present numerical 
predictions for the quality factors and mode volumes of 
localized defect modes as a function of geometric param- 
eters, and discuss some experimental challenges related 
to the coupling of a defect cavity to gas-phase atoms. 

1 Introduction 

A variety of passive and active optical devices can be con- 
structed by introducing point or line defects into a peri- 
odic array of holes perforating an optically thin semicon- 
ductor slab. In such structures, light is confined within 
the defect regions by the combined action of distributed 
Bragg reflection and internal reflection. This technique 
has been employed in making a microcavity semicon- 
ductor (InGaAsP) laser Q (emitting at A = 1.55 fim), 
and for demonstrating Si optical waveguides with sharp 
bends ||^ . One can thus easily envision the fabrication of 
integrated optical "networks" on a single chip, with nu- 
merous microcavity-based active devices linked by pas- 
sive waveguide interconnects. 

The combination of high quality factor and extremely 
low mode volume that should be obtainable in point- 
defect microcavities makes the photonic crystal (PC) 
paradigm especially attractive for experiments in cav- 
ity quantum electrodynamics (cavity QED) with 
potential applications in quantum information technol- 
ogy. The primary focus of this paper will be to dis- 
cuss the optimization of microcavity designs for cavity 
QED with strong coupling between defect modes and 
gas-phase neutral atoms. This new paradigm for cavity 
QED poses formidable technical challenges with regard 
to atom trapping and the characterization of surface ef- 
fects, so we include some discussion of these issues and 
of our current approaches to addressing them. 

In Section II of this article we begin by analyzing el- 
ementary microcavities formed by changing either the 
refractive index or radius of a single defect hole. We 
discuss in detail some problems in previous FDTD cal- 
culations done by our group as intuitions gained 



from the resolution of these problems are important in 
the microcavity design optimization for cavity QED. In 
Section III we present our designs of microcavities opti- 
mized for strong coupling between the cavity fleld and an 
atom trapped within a hole of the PC. Based on numer- 
ical analysis using the 3D flnite-difference time-domain 
(FDTD) method, we predict that quality factors over 
3 X 10'' can be achieved in these structures. We also 
discuss in this section the difficulty of estimating surface 
interactions between a trapped atom and the semicon- 
ductor substrate, and propose a technical strategy for 
confining gas-phase atoms within a defect microcavity. 
Finally, in Section IV we describe a procedure that we 
have developed for microcavity fabrication. 
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Figure 1 : Optically thin slab patterned with a hexagonal 
array of air holes. 



2 Single defect hole 

The microcavities analyzed in this paper are formed by 
introducing point defects into a dielectric slab patterned 
with a hexagonal array of air holes. The thickness of the 
slab is d and its refractive index is equal to 3.4. The 
spacing between holes is denoted by a and the hole ra- 
dius by r, as shown in Figure |l|. We use A to specify 
the optical wavelength in air. In our calculations of mi- 
crocavity quality factors, the boundary for separation of 
vertical from lateral loss (i.e. the vertical quality factor 
Q± from the lateral quality factor Q^) is positioned ap- 
proximately at A/2 from the surface of the membrane. 
As the number of PC layers around a defect is increased, 
Q|| increases and the total quality factor Q approaches 
Q±. We will adopt a coordinate convention in which 
X — 0, y = 0, z = denotes the center of the cavity and 
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z = is the middle plane of the slab. 

The simplest method of forming a microcavity within 
the structure shown in Figure]^ is to change the radius or 
index of refraction of a single hole. The former method 
is more interesting from the perspective of fabrication, 
since the lithographic tuning of geometric parameters of 
individual holes is a simple process, but in this section 
we will consider both methods. By increasing the ra- 
dius of a single hole an acceptor defect state is excited, 
i.e., pulled into the band-gap from the dielectric band. 
On the other hand, by decreasing the radius of an in- 
dividual hole (or by tuning its refractive index between 
1 and the refractive index of the slab) a donor defect 
state is pulled into the band-gap from the air band 
Acceptors tend to concentrate their electric field energy 
density in regions where the larger (semiconductor) re- 
fractive index was located in the unperturbed PC, while 
the electric field energy density of donors is concentrated 
in regions where there was air in the unperturbed PC. 
Since the electric field intensity in air regions is small in 
case of acceptor defects, these are not good candidates 
for strong coupling with a single gas-phase atom that 
would be trapped within a hole. In this article we will 
thus focus on donor states. For a discussion of acceptor 
states excited in an optically thin slab perforated with a 
hexagonal PC array, readers are referred to Ref. |6j. 
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Figure 2: Electric field intensity pattern of the x-dipole 
mode. Parameters of the structure are r/a = 0.3, d/a = 
0.6 and in the used discretization a = 15. The refractive 
index of the slab is 3.4 and of the defect is 2.4. The 
plotted intensity patterns are for the x — y plane at the 
middle of the slab. 



2.1 Changing the refractive index of a 
single hole 

Microcavity formation by alteration of the refractive in- 
dex of a single defect hole in a hexagonal PC has been 
analyzed previously by our group W. In that analysis. 



we predicted that dipole-like donor states (such as the 
x-dipole mode shown in Figure |^) with quality factors up 
to 3 X 10'' should exist. We now believe that the qual- 
ity factors of such microcavities are limited to several 
thousand, for reasons discussed below. This discussion 
reveals the extreme sensitivity of microcavity quality fac- 
tors to small distortions to the local PC geometry, which 
will later be used as a powerful design tool in optimiza- 
tion for cavity QED. 
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Figure 3: Deformation of the analyzed structure intro- 
duced by the application of mirror boundary conditions 
along the x and y axes and through the center of the 
defect. The applied mirror boundary conditions select 
the x-dipole mode, whose electric field pattern is shown 
in Figure ^ Holes on the x-axis (in the central row, con- 
taining the defect) are elongated by 1 point in the y direc- 
tion. The refractive index of the defect is rifi^ff.ct = 2.4. 
Parameters of the structure are r/a ~ 0.3, d/a ~ 0.6 
and in the used discretization a = 15. 



In our previous work, mirror boundary conditions were 
applied in the x, y and z directions to achieve an eight- 
fold reduction in the computational grid size. We have 
since realized that the manner in which even (symmet- 
ric) mirror boundary conditions are implemented in our 
finite-difference time domain (FDTD) code results in nu- 
merical output that properly corresponds to an analyzed 
structure with slight deformations relative to the in- 
tended structure. For example, the set of mirror bound- 
ary conditions used to select the x-dipole mode in a de- 
fect cavity leads to a deformation of the structure as 
shown in Figure || Holes on the x-axis (including the de- 
fect) are elongated in the y direction by 1 point but the 
distance between holes in x and y directions is preserved. 
Because hole-to-hole distances are preserved under this 
deformation, the half-spaces y > 1/2 and y < —1/2 actu- 
ally maintain the unperturbed PC geometry when holes 
in the central row are elongated by 1/2 points in both 
the ±2/ directions. The symmetry of the PC surround- 
ing the defect is therefore broken, and this contributes to 
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artificially high quality factors for x-dipole modes. An 
even mirror BC was also applied in the z direction in our 
previous analysis, causing a slight error in the thickness 
of the slab. The correct dja ratios of the structures an- 
alyzed in Ref. Q would be 0.6, 0.46 and 1, instead the 
values of 0.53, 0.4 and 0.93, as noted there. 

2.1.1 Perfectly symmetric cavity 

We re-analyzed the structure with r ja = 0.3, d/a = 0.6, 
ndefect = 2.4, the refractive index of the slab equal to 
n = 3.4, 5 layers of holes around the defect, and a = 15 
grid points. For this set of parameters, the predicted 
Qi_ in Ref. Q for the x-dipole mode was 3 x 10^. In the 
present analysis, even mirror boundary conditions were 
applied to the lower boundary in the z direction only, to 
reduce the computation size by one half and to eliminate 
TM-like modes. Absorbing boundary conditions were 
applied to all boundaries in x and y directions and to 
the upper boundary in z direction. Under these BC's the 
intended symmetry of the defect structure was achieved. 
The initial field distribution was chosen in such a way 
as to excite x or y dipole modes selectively. For the x- 
dipole mode we now calculate Q\\ — 2260, = 1730 
and ^ = 0.3137, and for the y-dipole mode we calculate 
Q|l = 1867, Qa_ = 1007 and f = 0.3182. 

The difference between parameters computed for the 
X and y dipole modes comes partly from the asymme- 
try of the structure introduced by imperfect discretiza- 
tion. In a 2D PC with infinite slab thickness, these two 
modes should be degenerate Q. In the thin slab how- 
ever, the y-dipole mode suffers more vertical scattering 
at the edges of holes and, therefore, has a lower Q±. For 
a non-symmetric applied initial field distribution leading 
to excitation of both x and y dipole modes, we calculated 
Q\\ = 2070, Qi_ = 1290 and f = 0.316. This mode can 
be represented as a superposition of the x and y dipole, 
with weighting factors depending on the initial field. Its 
quality factor and resonant frequency depend on these 
weighting factors. 
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Figure 4: (a) Q factors and (b) frequencies of x-dipole 
modes in the structure shown in Figure |^, as a function 
of the elongation parameter p. The structure parameters 
are r/a = 0.3, d/a = 0.6, n = 3.4, ndefect = 2.4, a = 15 
and 5 layers of holes surround the defect. 



2.1.2 Asymmetric cavity 

Our explanation for the discrepancy between our current 
and previous results leads us to ask whether a "real" (in- 
tentional) elongation of the central row of holes along the 
y-axis may actually improve the quality factor (Q) of the 
a;-dipole mode. A possible disadvantage of this approach 
is the excitation of additional states caused by enlarging 
holes in the central row. We have therefore analyzed how 
a/ A and Q of the x-dipole mode in this structure changes 
as a function of the elongation parameter p. Holes on the 
X-axis (including the defect) are elongated in both the 
±2/ directions by p/2 points, in such a way that hole-to- 
hole distances are preserved and the half-spaces y > p/2 
and y < —p/2 maintain the unperturbed PC geome- 
try. The structure parameters are r/a — 0.3, d/a — 0.6, 
n — 3.4, ndefect = 2.4 and 5 layers of holes surround the 



defect. The periodicity a used in these calculations is 
equal to 15 grid points, and the elongation step Ap = 1 
corresponds to a/15. The results are shown in Figure ^. 
It is interesting that the frequency of the mode decreases 
as p increases, even though the amount of low refrac- 
tive index material increases. However, the net amount 
of low refractive index material does not matter; what 
does matter is where the low refractive index is posi- 
tioned relative to the unperturbed PC. The explanation 
of the decrease in frequency is very simple, if we recall 
the x-dipole mode pattern shown in Figure ^. This is a 
donor type defect mode, which concentrates its electric 
field energy density in low refractive index regions of the 
unperturbed PC. As p increases, layers of PC holes are 
moved away from the defect in y direction. For example, 
the 71-th layer of holes parallel to the x-axis will be posi- 
tioned aX y — ±ria\/3/2 ±p/2, instead oi y — ±na \/3/2. 
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Figure 5: Q factor for p = 3, in the structure shown in 
Figure ^, as a function of number of PC periods around 
the defect. The structure parameters are r/a = 0.3, 
d/a = 0.6, n = 3.4, Udefect = 2.4 and a = 15. 
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Figure 6: Band diagram for TE-hke modes of a thin 
slab {n = 3.4, d/a — 0.75) surrounded by air on both 
sides and patterned with a hexagonal array of air holes 
(r/a = 0.275). 



Therefore, a large refractive index material will be po- 
sitioned at places where the mode expects to "see" air, 
leading to a decrease in mode's frequency. By tuning 
the mode's frequency across the band-gap, we can also 
tune its Q factor, as noted previously Q. For p = 3, 
Q± reaches the value of 20000. The analyzed x-dipole 
mode mostly resonates in the direction of y-axis, i.e. in 
the TX direction of photonic crystal. The tuning of the 
elongation parameter p is, therefore, analogous to tuning 
of a spacer in the micropost cavity, which leads to tun- 
ing of mode's frequency and Q factor. A more detailed 
explanation of the effect of elongations on Q factors is 
given in the appendix of this article. By increasing the 
number of PC periods around the defect, the total qual- 
ity factor Q approaches Q±, as shown in Figure |[ It is 
important to note that Q\\ does not increase exponen- 
tially with the number of PC layers around the defect. 
Instead, it saturates at large number of PC layers. The 
reason is in the choice of a boundary for separation of 
Q± from Q||, positioned at approximately A/ 2 from the 
surface of the membrane. From the radiation pattern 
of the x-dipole mode, we can see that some portion of 
the out-of-planc loss (mostly in the x direction) still gets 
collected in (3||. This loss cannot be suppressed by in- 
creasing the number of PC layers around the defect, and 
it determines the value at which Qn saturates. However, 
much larger out-of-plane loss is collected in Q±, which 
ultimately determines the total quality factor Q. 
Therefore, a dramatic improvement in Q factors of dipole 
modes over single defect microcavities can be obtained 
by introducing a novel type of PC lattice defect, consist- 
ing of the elongation of holes along the symmetry axes. 
We call it a fractional edge dislocation, by analogy with 
edge dislocations in solid state physics. Edge disloca- 
tions are formed by introducing extra atomic planes into 



the crystal lattice. On the other hand, we here insert 
only fractions of atomic planes along the symmetry axes 
of photonic crystal. 

Unfortunately, at this time, wc do not know how to 
control the refractive index of a single PC hole during 
the fabrication. For this reason, in the next subsection 
we will consider alternative methods of forming single 
defect microcavities, which are much easier to construct 
by microfabrication. 

Even after the elongation of holes on the x-axis by 1 
point, Q factors of 30000 were not obtained. This means 
that the application of mirror BC at a; = and y = 
planes causes additional effects that lead to the overes- 
timation of quality factors. One of the reasons may be 
that the excited dipole mode does not have a symmetry 
described by the applied discretized mirror BC's. This 
may be partly due to the structure imperfection caused 
by discretization. In order to avoid problems caused by 
BC's all the calculations in this paper were done by ap- 
plying the even mirror symmetry to z = plane only, in 
order to select TE-like modes and to reduce the compu- 
tation size by one half. Absorbing boundary conditions 
were applied to all boundaries in x and y directions and 
to the upper boundary in the z direction. To prove that 
the application of mirror BC at the lower z boundary 
does not change Q, we also analyzed entire structures 
with absorbing BC's applied to all boundaries and ob- 
tained the same results as in the analysis of one half of 
the structure. 

2.2 Reducing the radius of a single hole 

Microcavities analyzed in this subsection are formed by 
reducing the radius of a single hole to rdef ■ We calcu- 
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Table 1: Q factors of dipole modes excited in microcavi- 
ties formed by decreasing the radius of a single PC hole. 
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lated parameters of excited dipole modes for a range of 
microcavity parameters and results are shown in Table 
H For all tabulated results, 5 layers of holes surround 
the defect and a = 20. Calculated Q's are not very im- 
pressive, but they do provide us with a good starting 
point for further optimization. Why did we decide to 
use a relatively small r ja ratio for the PC? According 
to our calculations, increasing rja within the analyzed 
range leads to increases in the band gap and reduction of 
lateral losses. However, vertical scattering at the edges 
of holes also increases and drops. It is therefore im- 
portant to find an optimum rja which leads to small 
vertical losses but preserves good lateral confinement, in 
order not to increase the mode volume too much. 

The band diagram for TE-like modes of a thin slab 
(n = 3.4, dja = 0.75) surrounded by air on both 
sides and patterned with a hexagonal array of air holes 
[rja — 0.275) is shown in Figure ||. These PC parame- 
ters are used in most of the calculations in the next sec- 
tion. From the comparison of the dipole mode frequen- 
cies tabulated in Table |l| and the band diagram shown in 
Figure ||, it can be confirmed that as the cavity mode's 
frequency approaches the bottom of the air band, verti- 
cal losses decrease, but lateral losses increase, which in 
turn leads to an increase in the mode volume. Our goal 
in the next section will be to reduce vertical losses and 
improve Q factors even further, while preserving small 
mode volumes. For this purpose we will explore frac- 
tional edge dislocations. In Figure |[ one can observe 
that an increase in the elongation parameter p can be 
used to tune the factor of a mode, but also leads to 
a decrease in the dipole mode's frequency. This implies 
that by increasing p, the mode is pulled deeper into the 
band gap, away from the air band edge, which leads to 
its better lateral confinement . Therefore, we can simul- 
taneously achieve a reduction in vertical losses and an 
improvement in lateral confinement (i.e. an increase in 
both Qj^ and Q\\, and a reduction in the mode volume). 

3 Cavities for strong coupling 

In this section we consider the design of PC microcavities 
to achieve strong coupling between the cavity field and 
a single gas-phase atom, that is, an atom located in free 
space rather than contained as an impurity in the dielec- 
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Figure 7: Parameters of the x-dipole mode in a single 
defect structure {rja = 0.275, dja = 0.75, Vdefja = 0.2, 
n — 3.4 and a = 20) as a function of the elongation 
parameter p: (a) mo; (b) A^'o 



trie slab. Our long-term goal is to investigate photonic 
bandgap structures for single-atom cavity quantum elec- 
trodynamics in the strong coupling regime [|. For this 
purpose, the microcavity mode quality factor (Q) has 
to be as large as possible and the mode volume (Vmode) 
as small as possible. These two design rules are also 
followed when designing PC microcavities for semicon- 
ductor lasers. However, in a cavity for strong coupling, 
an atom must be trapped at the point where it interacts 
most strongly with the cavity field. Therefore, an ad- 
ditional design goal is imposed in this case: the cavity 
mode should have the E-field intensity as high as possi- 
ble in the air region. When designing a laser cavity, the 
problem is opposite: one tends to maximize the overlap 
between the gain region and the cavity field and, there- 
fore, wants to have the strongest E-field in the semicon- 
ductor region. 

Mode volume {Vmode), critical atom (A'o) and photon 



5 



(mo) numbers are defined as follows: 
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where n is the cavity field decay rate, proportional to 
the ratio of the angular frequency of the mode (wq) and 
the mode quality factor (Q): 
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7^ is the atomic dipole decay rate (2.6 MHz for Cesium) 
and g is the coupling parameter at the point where we 
want to put an atom: 
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go denotes the vacuum Rabi frequency: 
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Strong coupling is possible if both Nq and toq are 
smaller than 1. Therefore, in order to predict whether 
the strong coupling can occur, we must calculate upper 
limits of A^o and mo and compare them to 1. In other 
words, it is acceptable if calculated critical numbers are 
overestimated. As the number of PC layers around the 
defect increases, the total quality factor Q approaches 
Q±^ and Vmode drops due to the better lateral confine- 
ment. Hence, we can calculate iVo and mo by assuming 
Q = Q± and using Vmode calculated for 5 PC layers 
around the defect. 

For all calculations in this section, the refractive index 
of the slab is n = 3.4, 5 layers of holes surround the 
central hole and a = 20. The elongation step Ap — 1 
corresponds to a/20, i.e. 5% of the lattice periodicity 
a. The material and PC properties are chosen in such a 
way that cavities operate at A = 852 nm (the wavelength 
corresponding to the D2 atomic transition in ^^^Cs). 

3.1 Single defect with fractional edge 
dislocations 

Let us study microcavities formed by reducing the radius 
of a single hole and simultaneously applying a fractional 
edge dislocation (of order p) along the x-axis. It is impor- 
tant to note that all holes along the a;-axis are elongated 
in this process, including the defect hole of the reduced 
radius. We will calculate the dependence of the x-dipole 
mode properties on parameter p. Parameters of the un- 
perturbed PC are: r/a = 0.275, d/a = 0.75, a = 20 and 
the defect hole radius is rdef/a = 0.2. The critical atom 



(iVo) and photon (mo) numbers are shown as a function 
of the elongation parameter p for the x-dipole mode in 
Figure ^ Vmode does not change significantly with p and 
it is approximately equal to 0.1(A/2)'^ for all structures. 
When p increases, the frequency of the mode moves away 
from the band edge, towards the center of the band gap, 
Q|| increases and Q± peaks at the value of 1 x 10* for 
p = 2. 
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Figure 8: Electric field intensity pattern of the x-dipole 
mode excited in a microcavity formed by reducing the ra- 
dius of a single hole and simultaneously elongating holes 
on the X-axis by 2 points, (a) slice through the middle 
of the membrane, in the x — y plane; (b) x — z half-plane 
(z > 0). 



From the electric field intensity pattern of the x-dipole 
mode shown in Figure ||, one can see that the elec- 
tric field intensity is very strong within the defect hole. 
Therefore, an atom trapped there should interact very 
strongly with the cavity field. From the calculated criti- 
cal atom and photon numbers, it then should be possible 
to achieve very strong coupling. At A = 8527T.m, the pa- 
rameters of such a cavity are r = 70nm, d = IQOnm, 
a — 250 and rdef = 50nm. Due to extremely small 
mode volumes in these cavities, strong coupling is possi- 
ble even for moderate values of Q, as Q± did not exceed 
^ 1 X 10* in the parameter range of Fig. 0. Further- 
more, mo is much smaller than Nq, which means that 
we can try to improve Q factors further at the expense 
of increasing Vmode- Q factors above 1 x 10* and sim- 
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ilar values of tuq and A'o can also be obtained for the 
cavity consisting of a single defect with Vdef/a = 0.2 
and a fractional edge dislocation of order p, produced 
in the photonic crystal whose parameters are r/a = 0.3, 
d/a — 0.65 and n — 3.4. 




Figure 9: Tuning 4 holes closest to the defect in TJ 
directions. Their radii are reduced to ri and they are 
simultaneously moved away from the defect in TJ direc- 
tions by r — ri . The radius of the central hole is r2 . 



3.2 Tuning holes around the defect 

Dipole modes are particularly sensitive to the geometry 
of holes closest to the defect. By tuning these holes, we 
can induce frequency splitting of dipole modes and dra- 
matically influence their Q factors. In Rcf. the vari- 
ation of two nearest neighbor holes along the x-axis was 
analyzed. Here, we will test the influence of changing 4 
holes closest to the defect in TJ directions. The analyzed 
structure is shown in Figure ^ The radius of the central 
hole is reduced to r2 and the radii of the 4 closest holes 
in TJ directions are reduced to ri. These 4 holes are 
simultaneously moved away from defect, by r — ri in TJ 
directions, which preserves the distance between them 
and the next nearest neighbors in the same directions. 
This design will improve the Q factor of the y-dipole 
mode and spoil the Q of the a;-dipole mode. We analyzed 
structures with various parameters, but our best result 
was obtained for r/a = 0.275, d/a = 0.75, r2/a = 0.2, 
Ti/a = 0.225 and a = 20. The electric field intensity 
pattern of the excited y-dipole mode is shown in Figure 
and its calculated parameters are a/A = 0.289 and 
= 4890. From Table we can see that a dipole mode 
excited in a single defect microcavity with this a, r/a, 
d/a has = 2078. Therefore, the tuning of 4 holes can 
lead to a substantial increase in Q of the x-dipole mode. 
The disadvantages of this design include the excitation of 



defect modes other than dipoles (coming from variation 
of several holes). 

Let us try to improve Q of this cavity even further, by 
also employing the idea of elongation of holes along de- 
sired directions. Our mode of choice is the y-dipole and 
we will elongate holes sitting on the y-axis by p points 
in the x direction in such a way that the half-spaces 
X > p/2 and x < —p/2 remain the unperturbed PC ge- 
ometry. The dependence of A'o (which decreases with 
Q) on parameter p is shown in Fig. V„iode does not 
change significantly with p and is in the range between 
0.09(A/2)3 and 0.12(A/2)3. The calculated toq is around 
5 • 10~^ for all structures. Again, a very strong coupling 
is achievable by this design, and we note that Q±_ at the 
point p = 2 achieves a value ~ 3.3 x 10**. 
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Figure 10: Electric field intensity pattern of the y-dipole 
mode excited in the cavity where 4 holes closest to the 
defect in TJ directions are tuned. Their radii are reduced 
to ri/a = 0.225 and they are simultaneously moved away 
from the defect in TJ directions by r — ri . The radius 
of the central hole is r2/a = 0.2. PC parameters are 
r/a = 0.275, d/a = 0.75, a = 20 and n = 3.4. 



3.3 Atomic physics 

Given that microcavities with strong coupling parame- 
ters can be designed and fabricated, two further tech- 
nical issues must be addressed in order to establish the 
feasibility of cavity QED with neutral atoms in PC's. 
First, we must identify a method for stably trapping an 
individual atom within a hole of the PC. Second, ac- 
curate estimates or measurements must be made of the 
surface interaction between such a trapped atom and the 
semiconductor substrate. Although definitive solutions 
to these challenges are the subject of future work, we in- 
clude in this subsection a brief discussion of each topic. 
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In addition to the creation of photonic crystals, mod- 
ern micro-fabrication techniques enable the patterning 
of either conductive wires or ferromagnetic materials at 
the micron scale and below. As a result, it should be 
possible to construct magnetic microtraps with field cur- 
vatures - 10^ G/cm2 § II 0. With this magnitude 
of field curvature an loffe trap could hold a Cs atom with 
Lamb-Dicke parameter rj ^ 0.035 in the radial direction, 
yielding Ax ~ 10 nm in the ground state of the trapping 
potential and enabling resolved-sideband laser cooling as 
a means of putting single atoms in the ground state. 




''o 0.5 1 1.5 2 2.5 3 3.5 4 

P 

Figure 11: Nq of the a;-dipole mode in the structure 
where 4 holes in TJ directions are tuned (as shown in 
Figures | and |o|), as a function of the elongation pa- 
rameter p. Holes on the y-axis are elongated by p points 
in the x direction in such a way that the half-spaces 
X > p/2 and x < — p/2 remain the unperturbed PC ge- 
ometry. 

We are currently investigating fabrication and laser 
cooling techniques for an atom-trapping scheme in which 
micron-scale wires would be deposited on the surface of 
the PC semiconductor substrate, such that the circular 
wire pattern of an loffe microtrap are arranged concen- 
trically around a defect microcavity. The trap designs 
discussed in Q have sufficiently large inner diameter not 
to disturb the photonic bandgap structure of the defect 
cavity. Such a wire arrangement would project a mag- 
netic field with a stable minimum at the geometric center 
of the microcavity, such that one or more atoms could 
be confined within the defect hole and would therefore 
experience strong coupling. A similar geometry could be 
envisioned for microtraps based on permanent magnets 
rather than current-carrying wires, which would have 
significant advantages in terms of heat load to the PC 
substrate. 

A neutral atom trapped within a hole of a PC struc- 
ture will experience surface (van der Waals) interactions 
that are quite difficult to estimate. The significance of 
such interactions is twofold. The sensitive dependence 



of the ground state energy shift on an atom's position 
will lead to mechanical forces that must be compensated 
by the trap design. Differential shifting of the atomic 
ground and excited states on the cavity QED transi- 
tion must also be accounted for, as this will introduce 
a position-dependent detuning relative to the fixed mi- 
crocavity resonance. Although van der Waals shifts can 
be computed for alkali atoms near dielectric or metal- 
lic boundaries with simple symmetry |TT[ [l^ , the case 
of an atom in a PC hole is far more complex. The lo- 
cal geometry seen by a trapped atom will be that of a 
cylindrical hole with finite extent, and a proper calcula- 
tion must take into account the overall modification of 
vacuum modes due to the extended photonic crystal. 



Figure 12: Electric field intensity patterns of the cou- 
pled dipole modes: (a) constructively; (b) destructively 
coupled defect states. 
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Furthermore, the resonant frequencies of many atomic 
transitions that connect to low-lying states, and there- 
fore contribute strongly to their van der Waals shifts, are 
above the bandgap of the semiconductor substrate. Such 
transitions will "see" an absorptive surface while those 
below the bandgap will see a dielectric surface. This set 
of factors brings the complexity of the desired calcula- 
tion well beyond that of existing analytic results in the 
literature. It should be noted that experimental mea- 
surements of the van der Waals shifts in our proposed 
system would be of significant interest for the general 
subject of quantum electrodynamics of semiconductors. 

We are pursuing a numerical strategy for estimating 
the magnitude of surface interactions. In a linear re- 
sponse approximation, it should be possible to compute 
leading-order contributions to the van der Waals shifts 
from FDTD simulations of the electromagnetic field cre- 
ated by an oscillating dipole source in the photonic crys- 
tal structure. While the nature of the code does not 
allow us to compute directly the backaction of the scat- 
tered field on the source dipole, we believe that elemen- 
tary field theory can be used to relate the simulated field 
to van der Waals shifts. Our findings will be reported in 
a forthcoming publication. 
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Figure 13: Electric field intensity patterns of the con- 
structively coupled dipole modes in the structure with 
the following parameters: rdef/a = 0.2, d/a = 0.75, 
r/a = 0.275, n = 3.4 and a = 20. Holes in columns 
containing defects are elongated by 2 points in the x 
direction 



3.4 Coupled dipole defect modes 

The significance of surface effects that could perturb 
atomic radiative structure within the small defect hole 
is still unknown. For that reason, wc will try to investi- 



gate ways of increasing the radius of the hole where the 
coupling between the atom and the cavity field should 
occur. Let us now analyze the cavity design where the 
strong E-field intensity can be achieved in the center of 
an unperturbed hole. The idea is to use coupling of two 
dipole defect states. Resonant modes of the microcavity 
formed by coupling two single defects are presented in 
Figure |l2|. Based on the resultant electric field intensity 
in the central, unperturbed hole, we call them construc- 
tively or destructively coupled defect states. They have 
different frequencies, as well as Q factors. We will ana- 
lyze properties of the constructively coupled state, since 
the central, unperturbed hole would be a good place for 
an atom. 

o o o o o 
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Figure 14: Electric field intensity patterns of the con- 
structively coupled dipole modes in the structure with 
the following parameters: rdef/a = 0.2, d/a — 0.75, 
r/a = 0.275, n = 3.4 and a = 20. Holes in rows contain- 
ing defects are elongated by 2 points in the y direction. 

We analyzed a series of structures with different pa- 
rameters. The best results were obtained for two cou- 
pled defects with r^ef/a = 0.2 in a PC with the fol- 
lowing parameters: d/a = 0.75, r/a — 0.275, n ~ 3.4 
and a = 20. Holes in the FX direction, in columns 
containing defects, are elongated by 2 points in the x di- 
rection. The mode pattern of the constructively coupled 
defect state is shown in Figure Parameters of the 
mode are: a/A = 0.29, Q± = 6100, Vmode = 0.19(A/2)3, 
mo = 1.5 • 10^^ and iVo — 0.0135. An atom can now 
be trapped in the central hole of unperturbed radius. 
For A = 852nm, this radius is r = 68nm, which is a 
significant improvement over the previous design, where 
an atom must be trapped within a 50nm radius hole. 
Again, a strong coupling is achievable in this cavity. 

An alternative way of forming the coupled defects 
state is represented in Figure |lj. We use the same PC 
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parameters as previously: rdf,f/a — 0.2, d/a — 0.75, 
r/a — 0.275, n = 3.4 and a = 20. Holes in rows contain- 
ing defects, are elongated by 2 points in the y direction. 






Figure 15: SEM pictures showing the top views of the 
fabricated structures. 



4 Fabrication 

We have recently developed the fabrication procedure 
for making these cavities in AlGaAs. The material and 
PC properties are chosen in such a way that cavities 
operate at A = 852nrn (the wavelength corresponding to 
the atomic D2 transition in ^^^Cs). 

The fabrication process starts by spinning of lOOrtm 
thick high molecular weight PMMA (polymethyl- 
methacrylate) on top of the wafer. The PMMA layer is 
subsequently baked on a hot plate at 150°C for 20 min- 
utes. A desired 2D PC pattern is beamwritten on the 
PMMA by electron beam lithography in a Hitachi S-4500 
electron microscope. The exposed PMMA is developed 
in a 3:7 solution of 2-ethoxyethanol:methanol for 30 sec- 
onds. The pattern is then transferred into the AlGaAs 
layer using the CI2 assisted ion beam etching. After 
that, the sacrificial AlAs layer is dissolved in hydroflu- 
oric acid (HF) diluted in water. HF attacks AlAs very 
selectively over Al^Gai-xAs for x < 0.4 jl^. There- 
fore, the percentage of Al in our AlGaAs layer is around 
30%. Finally, the remaining PMMA may be dissolved in 
acetone. 

Three SEM pictures showing top views of fabricated 
microcavity structures are shown in Figure |l^. We are 
currently working to measure the passive optical prop- 
erties of such microcavities in order to validate our the- 
oretical predictions. 

5 Conclusion 

In conclusion, we have theoretically demonstrated that 
PC cavities can be designed for strong interaction with 
atoms trapped in one of PC holes. At present, we are 
working on further optimization of the design and the 
characterization of fabricated structures. 

Critical issues for further investigation include efficient 
coupling of light in and out of the PC microcavity, as 
well as accurate estimation of surface effects that could 
perturb atomic radiative structure within the small de- 
fect hole. The extremely small mode volume in these 
structures also poses an interesting theoretical question 
of how standard cavity QED models must be modi- 
fied when the single-photon Rabi frequency exceeds the 
atomic hyperfine spacing. 



The mode pattern of the constructively coupled defect 
state is shown in Figure 14. Parameters of the mode 
are: a/A = 0.288, Qj_ = 12120, V^ode = 0.14(A/2)3, 
mo = 1.4 • 10^'' and A^o = 0.0063. Strong couphng is 
achievable for an atom trapped in any of the two central 
holes of the unperturbed radius (positioned between the 
defects). For A = 852nm, this radius is again r ~ 68nm. 



6 Appendix: The effect of frac- 
tional edge dislocations 

Any wavefront can be considered as a source of secondary 
waves that add to produce distant wavefronts, accord- 
ing to Huygens principle. Let us assume that we know 
the field distribution accross the plane S, positioned in 
the near field, above the free standing membrane and in 
parallel to the membrane surface. The far fields can be 
considered as arising from the equivalent current sheets 
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on this plane. Therefore, we can calculate the far field 
distribution and the total averaged radiated power into 
the half-space above the plane 5'|l4|: 









2 




+ 




V 







(8) 




where N and L represent radiation vectors, whose 
components in the rectangular coordinate system are 
proportional to Fourier transforms of tangential field 
components at the plane 5*1^51: 



= -FT2{Hy) 
Ny - FT2{H,) 
Lx = FT2{Ey) 



-FT2{Ex) 



27r 



FT2{f(x,v)) 



— sinOyxcosq) + ysinq 
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(10) 

(11) 
(12) 

(13) 
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the improvement in the Q factor after the application of 
fractional edge dislocations is due to the suppression of 
the wave- vector components composing the defect mode, 
which are located within the light cone. A more detailed 
explanation of this phenomenon and how it can be used 
to improve Q factors of other types of modes will be pre- 
sented in our forthcoming publication [p^ . 

The Q factor of a mode can be expressed as Q = , 
where W is the total energy of a mode in the half-space 
z > 0. The comparison between the Q factor calculated 
using the method presented in this appendix, and (5± 
previously estimated using the FDTD, for the structure 
from the subsection 2.1.2, is shown in Figure |l6|. The 
plane S is positioned directly above the surface of the 
membrane in this case. From Figure |^ it follows that 
the total Q factor saturates at about 17000, when the 
number of PC layers around the defect increases in the 
structure with p = 3. This is very close to the maximum 
Q value estimated from the expression ^. 
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Therefore, just by knowing the Fourier transforms of 
the tangential field components at the plane S*, we can 
evaluate the time averaged radiated power P. From the 
previous expressions it is clear that the wave-vector of 
interest fc|| lies within the light cone for any values of 
angles 9 and in the circular polar coordinate system 
(i.e. |fc||| < ^). This implies that the radiated power 
P depends only on the wave-vector components located 
within the light cone. It is also clear that by suppressing 
the Fourier components within the light cone, one can 
reduce P. Having in mind that for the x-dipole mode Ey 
and Hx fields are odd with respect to both x and y sym- 
metry axes, Ny and Lx do not contribute significantly to 
the integral in expression || (they are both equal to zero 
at any point in the k-space with either kx or ky equal to 
zero). On the other hand. Ex and Hy field components 
are even with respect to both x and y axes and their 
Fourier transforms are generally non-zero at small wave- 
vector values. However, by tuning the elongation factor 
p, one can balance the energy in the positive and nega- 
tive field lobes and minimize the Fourier components of 
Ex and By within the light cone. This also leads to a 
decrease in the radiated power P. We can conclude that 



Figure 16: Q factor computed using the FDTD method, 
and from the expression |^, for the structure analyzed in 
the subsection 2.1.2. 
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